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Final State Interactions effects are discussed in the context of Monte Carlo simulations of neutrino- 
nucleus interactions. A role of Formation Time is explained and several models describing this effect 
are compared. Various observables which are sensitive to FSI effects are reviewed including pion- 
nucleus interaction and hadron yields in backward hemisphere. NuWro Monte Carlo neutrino event 
generator is described and its ability to understand neutral current tt" production data in ~ 1 GeV 
neutrino flux experiments is demonstrated. 



I. INTRODUCTION 

New generation of neutrino oscillation parameters 
measurements require a good knowledge of neutrino- 
nucleus cross sections. Experimental data analysis is al- 
ways based on predictions from Monte Carlo (MC) event 
generators JH ■ In the 1 GeV energy region, characteristic 
for several oscillation experiments (MINOS, T2K, Mini- 
BooNE, NOvA) the use of the Impulse Approximation 
(lA) picture in which neutrinos scatter on individual 
quasi- free nucleons is well justified. In this picture any 
neutrino-nucleus interaction becomes a two-step process: 
(i) the primary scattering on a bound nucleon, and (ii) 
Final State Interactions (FSI) affecting the hadrons pro- 
duced at the step (i). The FSI contribute significantly 
to the systematic errors in neutrino oscillation measure- 
ments so it is important to develop models to describe 
them better and also to understand the models' limita- 
tions fs!]. 

MC codes used in major neutrino oscillation experi- 
ments (FLUKA (3, NUANCE [|, NEUT Q, GENIE 
) in their description of FSI effects rely on the model 
of intra-nuclear cascade (INC) [1]. It is a semi-classical 
approach in which some quantum effects can also be in- 
corporated (Pauli blocking, formation time (FT) , nucleon 
correlations). Theoretical arguments for the applicability 
of the cascade model go back to the works of Glauber 9] . 
More recently the investigation of the cascade model in 
the A resonance region was done in ■ The model pre- 
dictions agree with the experimental data for the pion- 
nucleus reaction cross sections, including the pion ab- 
sorption. 

While the basic idea behind the models of FSI in the 
MC codes is always the same, numerical implementations 
are quite different reflecting priorities of particular neu- 
trino experiments (target, detection technique etc). 
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An important and not sufficiently understood ingredi- 
ent in the INC models are Formation Time (FT) effects. 
On the most fundamental level the FT is related to the 
Quantum Chromo dynamics phenomenon called the color 
transparency (CT), proposed by Brodsky and Mueller 
For high enough four-momentum transfers a quark 
system is created with a small transverse size (point- 
like configuration - PLC) which is supposed to suppress 
hadrons re-interactions. As the typical size of the PLC is 
of the order of 1/|(5| fl^], the CT effects are expected to 
be seen mostly at higher energies. Moreover, two-quark 
systems are more likely to create PLC than three-quark 
ones so the effect is expected to be larger for pious, than 
for nucleons. 

Independent phenomenological considerations [13|, 
[3] led to the construction of approximate models of 
FT. As will be shown in Sect IIVI many evaluations of 
basic parameters which determine the size of FT effects 
have been proposed. It seems important to study them 
explicitly in the context of neutrino measurements. For 
example, the FT effects are in the obvious interplay with 
the pion absorption, reducing its probability in a non 
trivial momentum dependent way. 

A validation of FSI models can be based on any 
hadronic observables as all of them are FSI sensitive. 
Such observables include: distributions of numbers of re- 
constructed hadron tracks, spectra of hadrons in the final 
state, their angular distributions etc. FSI models used in 
neutrino MC simulations can also be validated on electro- 
and photo- nucleus observables. In the analysis of the 
NOMAD high energy neutrino scattering data J^] the 
introduction of FT was necessary to get the agreement 
with the measurements of backward moving protons and 
pious. It is interesting that the analysis of hadron- nucleus 
scattering data within INC models indicate that also at 
lower energies the FT effects can be quite important fl^ . 

In this paper FSI effects are modeled within the NuWro 
MC event generator 17]. NuWro covers neutrino energy 
range from a few hundreds MeV (the limit of applicability 
of the Impulse Approximation - lA) to several TeV. The 
code has flexibility to include Spectral Function [l^ for- 
malism with sophisticated nuclear effects, as an alterna- 
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FIG. 1. A block diagram of NuWro INC algorithm. 



tive to the Fermi Gas model or an momentum dependent 
effective potential 19]. NuWro allows for comparisons 
to the data reported by experimental groups in the FSI 
effects included format. 

We will consider a model of FT which is validated on 
the NOMAD backward moving pions data. We will then 
discuss the NC 7r° production data and see how impor- 
tant the FT effects are for the understanding of the ex- 
perimental data. NC lir^ is a very important process 
because it is a background to — )> oscillation search 
in water Cherenkov detectors: it can happen that one of 
two photons from 7r° decay remains undetected and the 
other is reconstructed as an electron. NC Itt" production 
is also a very useful reaction to validate FSI models in the 
1 GeV energy region. It is very sensitive to pion absorp- 
tion and it is important to investigate how relevant are 
FT effects making the nuclear environment more trans- 
parent for pions produced inside nucleus. 

The paper is organized in the following way. In Sect. 2 
a general description of the NuWro MC model is given. 
In Sect. 3 the NuWro FSI model based on the theoreti- 
cal approach of Oset [l^l is described. Several tests are 
reported showing a good agreement with the original nu- 
merical implementation. Sect. 4 contains a sumary of 
various ways to model the FT. Various approaches con- 
sidered in the context of neutrino interactions and pa- 
rameters used in theoretical computations and in MC 
codes are discussed. In Sect 5 the NuWro predictions are 
compared with the NC Ivr" production data and the sig- 



nificance of the FT effects is discussed. Our conclusions 
are contained in Sect 6. 



II. NUWRO 

NuWro is a neutrino event generation software devel- 
oped at the Wroclaw University 17[. The main motiva- 
tion for the NuWro authors was to have a tool to investi- 
gate the impact of nuclear effects on directly observable 
quantities with all the FSI effects included. Since 2005 
it evolved into a fairly complete neutrino interactions 
modehng tool. Its basic architecture is similar to better 
known MCs like NEUT or GENIE. All major neutrino- 
nucleus interaction channels are implemented and the 
commonly used relativistic Fermi Gas (EG) model is for 
certain nuclei replaced with the more realistic Spectral 
Function model [18]. The NuWro FSI code has recently 
been updated by implementing the Oset model of ef- 
fective pion-nucleon cross sections and several options for 
the FT. Other upgrades include: parameterization of the 
multipion production cross section in pion-nucleon colli- 
sions based on the available data and the implementation 
of angular distributions in elastic and charge exchange 
pion-nucleon scattering based on the SAID model [20] . 

With the inclusion of realistic beam models and a 
detector geometry module NuWro is becoming a fully- 
fledged MC event generator ready for use in neutrino 
experiments. 



A. Interactions 

In NuWro there are four basic dynamic channels: 
quasi-elastic (QEL), resonance (RES), more inelastic 
(DIS) and coherent pion production (COH), each can be 
either in the charged current (CC) or in the neutral cur- 
rent (NC) mode. The eight channel/mode combinations 
can be individually enabled or disabled. 

For each (but the coherent) channel a particular nu- 
cleon which will take part in the interactions is picked up 
with nuclear matter density used as the probability den- 
sity. Its momentum is chosen from a ball with the radius 
set to the Fermi momentum (or the local Fermi momen- 
tum calculated for that density in the case of LDA) or 
obtained as a draw from the Spectral Function. 



1. QEL 

The CC quasi-elastic and NC elastic reactions are han- 
dled by the QEL channel. It uses the standard Llewellyn 
Smith formulae [21] with several options for the vector 
form factors (dipole, BBA03 [22], BBBA05 Alberico 
et al HI). 

The global and local relativistic EG models or SF ap- 
proach are typically used but the kinematics based on the 
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FIG. 2. Probability (per fermi) of microscopic pion-nucleon 
interactions as a function of a distance from the centre of an 
iron nucleus. Pion kinetic energy Tj; = 165 MeV. The Oset 
model results are taken from [id ]. 



Axial form-factors are taken from the reanalysis of the 
ANL and BNL bubble chambers data [Hi. Not using 
the standard resonance production Rcin-Sehgal model 
[29 II to describe a contribution from higher resonances 
is justified by the quark-hadron duality hypothesis (soj 
and by the fact that higher resonances cannot be sep- 
arated in lepton-nucleus scattering. The non-resonant 
background is modeled as a fraction of the DIS contribu- 
tion for W G (1.3, 1.6) GeV scaled so as the passage to 
the pure DIS channel be smooth. 



3. DIS 

The DIS channel is defined as > 1.6 GeV. The to- 
tal cross sections are evaluated using the Bodek-Yang 
prescription [3l| . The Pythia6 hadronization routine 
is called for specific quark configurations [3^ to allow 
their meaningful use also in the small W region down to 
1.2 GeV. 
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FIG. 3. Probability of macroscopic quasi-elastic or absorp- 
tion interactions as a function of an impact parameter b for 
tt"'""' Ca scattering with pion kinetic e nerg y Tk — 180 MeV. 
The Oset model results are taken from [lOll. 



momentum dependent nuclear potential [l^ is also avail- 
able. Currenly, the Spectral Functions for carbon, oxy- 
gen, argon, calcium and iron are implemented in NuWro 
according to the tables obtained from Omar Benhar or 
calculated in [25|. In the SF mode the de Forest kine- 
matical prescription [2^ is used. 



2. RES 

The RES channel is defined as < 1.6 GeV, where 
W is the invariant hadronic mass. The dominant contri- 
bution comes from the single pion production mediated 
by the A(1232) resonance according to the model [27| . 



4. COH 

The coherent pion production is implemented using the 
Rein-Sehgal model [33| with lepton mass corrections. 

III. NUWRO FSI MODEL 

The NuWro FSI effects are described in a framework 
of the INC model Q. The neutrino interaction point is 
selected inside nucleus according to the nuclear matter 
density. All secondary hadrons propagate through nu- 
cleus and can interact with nucleons inside. In the code 
the 0.2 fm step length is assumed. For smaller values of 
the step the results remain the same, only running time 
increases. Between the collisions hadrons are assumed to 
be on-shell and move in straight lines. At each point of 
their path it is decided if there was an interaction or not. 
This is done based on an effective cross section model. 
The generic re-interaction algorithm is independent on 
the dynamics used, also different models of nuclear den- 
sity can be used. The particular dynamics is taken from 
the Oset model [k?!, as it has solid theoretical founda- 
tions. The model is supposed to work well in the most 
important A region for pion kinetic energies in the range 
85 — 350 MeV. Outside this region the cross sections are 
obtained from parameterizations of the available pion- 
nucleon cross section data. 

The basic FSI scheme consists in putting nucleons and 
pions produced in the primary and also in secondary in- 
teractions to a queue and repeating the following until 
the queue gets empty: 

(a) take a particle from the queue, 

(b) examine the nucleus density at its position. 
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FIG. 4. TT C absorption cross section. The data points are 
taken from: Ashery [s^, Navon [3^ . Jones [13] and Giannelli 
[3^ . The solid line shows NuWro predictions. 



FIG. 6. tt'*' ^^C charge exchange cross section. The data 
points are taken from: Ashery [ssl ]. Navon [s^, and Jones 
[33 • The solid line shows NuWro predictions. 
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FIG. 5. tt"*" ^^C inelastic cross section. The data points are 
taken from: Ashery 35], Jones [s^, and Levenson [39|. The 
solid line shows NuWro predictions. 



(c) calculate the mean free path, 

(d) probe the exponential distribution for the particle 
paths, 

(e) if the selected path is bigger then 0.2 fm adjust the 
particle position by 0.2 fm, 

• if the particle is still in the nucleus put it at the 
end of the queue, 

• if the particle is outside the nucleus put it to the 
list of outgoing particles, 

1° nucleons kinetic energy is diminished by the 
value of the potential: 



V = Jm2 



M ■ 



MeV 



where kj is the Fermi momentum and its 
momentum adjusted so that it remains on 
the mass shell. 
2° if nucleons kinetic energy is smaller then V 
the step 1° cannot be completed. The nu- 
cleon is assumed to be unable to leave nu- 
cleus. It is reinserted to the nuclear matter 
and its kinetic energy contributes to the nu- 
cleus excitation energy, 

(f) if the selected path is smaller that 0.2 fm assume that 
interaction with nuclear matter happened at this very 
place 

(g) probe the target nucleon momentum from the Fermi 
ball with the local Fermi momentum calculated from 
the density at that point 

(h) select the type of interaction and generate the kine- 
matics 

(i) check if none of the resulting nucleons is Pauli 
blocked; in the case of Pauli blocking forget the in- 
teraction, reinsert particle to the queue at the failed 
interaction point (Pauli blocking effects can be also 
included by means of increased values of the mean 
free paths and then this step of the algorithm must 
be skipped) 

(j) if the interaction was not Pauli blocked, all the parti- 
cles in the final state are put to the end of the queue. 

(k) if FT/FZ effects apply also to secondary interactions 
(this is model dependent) the particles positions are 
accordingly adjusted. 

The nucleus radius is defined as a distance from the 
center, where the density is smaller by a factor of 10^ 
than the maximal one. 
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As a result of some nucleons joining the INC, the nu- 
clear matter density is reduced but the shape of the den- 
sity profile is assumed to be unchanged. 

A. The Oset model 

On the microscopic level the Oset model includes the 
quasi-elastic pion-nucleon reaction (including the charge 
exchange channel) and the pion absorption with two- 
and three-body absorption mechanisms. The interaction 
probability per time unit is: 



Pdt = lmCn)dt = -2 lm(Vopt)dt (1) 

where w is the pion energy, 11 is the pion self-energy, 
and Vopt is the optical potential. 

In the simplest case of n'^p — > ir'^p p-wave scattering 
calculations lead to the result: 

where /* is ttNA coupling constant (/*^/47r ~ 0.36), 
is the pion mass, qc.m. is the pion momentum in the 
centre of mass system, Ga is A propagator, F its width 
and pp is proton density. 

An important in-medium effect is the A self-energy. 
Its imaginary part can be parameterized as [s^ ]: 

ImI]A(c^) - - [Cgip/por + CA2{p/pof + CA3{p/poy] 

(3) 

The A width is modified ^ ^ ImEA, changing 
the A propagator and producing extra terms in Eq. ([2|). 
proportional to functions C's present in Eq.(|31). The term 
proportional to Cq corresponds to higher order quasi- 
elastic scattering and the terms with Ca2 and Ca3 corre- 
spond to two- and three-body absorption, p is the nuclear 
matter density and po = 0.17 fm~^ is the normal density. 

The final expression for the interaction probability in 
the nuclear matter is: 

P = - f yf^^wlf— )^c.™.|GA(<7+fc)pJf(,+fc) 

uj J [2Try 3 J 2 

(4) 

where n{k) is the occupation number for pro- 
tons/neutrons. 

The A self-energy depends strongly on the nuclear den- 
sity and the pion absorption is more likely to occur in the 
central part of the nucleus. 

Finally, one introduces improvements to the model 
coming from: the ttN interaction s-wave contribution, 
the real part of the optical potential and finite size ef- 
fects. 
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n=l n—2 n—3 
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0.90 0.09 0.01 
0.89 0.10 0.01 


0.69 0.25 0.05 0.01 
0.67 0.24 0.07 0.02 



TABLE I. Probabilities that macroscopic quasi-elastic process 
proceeds through n microscopic collisions. Oset model results 
are taken from p^ . 
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0.81 0.17 0.02 
0.87 0.12 0.01 


0.37 0.41 0.17 0.04 
0.41 0.37 0.16 0.05 



TABLE II. Probabilities that pion absorption occurs after nth 
quasi-elastic microscopic scatterings. Oset model results are 
taken from 



B. NuWro implementation of the Oset model 

The Oset model is implemented in NuWro by means 
of tables containing the cross-sections as the functions 
of pion kinetic energy at various nuclear matter densi- 
ties. Because finite size effects are not universal it was 
necessary to prepare tables for each isotope separately. 

In the analysis of the performance of the cascade model 
one should distinguish microscopic (pion-nucleon) and 
macroscopic (pion-nucleus) reactions. Figs [2][3] show a 
comparison between our implementation and the original 
Oset model. Fig. [5] shows the inverse of mean free paths 
(or equivalently: interaction probabilities per fermi) for 
two microscopic interactions as functions of the distance 
from a nucleus center. In both cases contributions from 
pion-proton and pion-neutron reactions are added. There 
is a significant dependence on the nuclear density: in 
the nucleus central region the absorption probability is 
large. In the peripheral region the quasi-elastic scattering 
dominates overwhelmingly. Fig. [3] shows probabilities of 
macroscopic processes as functions of the impact param- 
eter. One can see that for small values of the impact 
parameter the absorption is more likely than the quasi- 
elastic scattering. One should remember that an incident 
pion can be absorbed also after one or more microscopic 
quasi-elastic scatterings. 

Tables Ulllll help to understand other aspects of 7r+ ^^C 
scattering. Again, we compare results from the origi- 
nal Oset paper with the NuWro implementation. Table U 
presents the probabilities that a macroscopic quasi-elastic 
event proceeds through n collisions. Table |ll] contains the 
probabilities that absorption occurs exactly after n mi- 
croscopic quasi-elastic pion scatterings. The results are 
shown for two values of incident pion kinetic energy: 85 
and 245 MeV. More energetic pious are likely to undergo 
several scatterings in which they loose a fraction of their 
energy until they fall into the absorption peak in the A 
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region. 

With a satisfactory agreement for microscopic ingre- 
dients of the Oset model, we present a comparison of 
the NuWro cascade model predictions with experimental 
data for the 7r+ ^^C scattering (Figs H- El). 

One can distinguish the following macroscopic pion- 
nucleus reactions: elastic, charge exchange, absorption, 
and inelastic scattering. We do not show the results for 
the double charge exchange reaction because the cross 
section is very small. For larger energies the inelastic 
cross section contains a pion production component. 

The cross section measurement of the charge exchange 
and absorption processes is straightforward. The inelas- 
tic cross section is obtained in the indirect way as: 
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^inel — ^total '^elastic - {c absorption + CTcEx), (5) 

where the elastic pion-nucleus cross section contribu- 
tion is evaluated based on theoretical and experimental 
arguments (for the details see [s^). 

The NuWro predictions are obtained in the standard 
way by arranging a homogeneous flux of pions and count- 
ing the particles in the final state assuming that at least 
one microscopic interaction took place. In the simu- 
lations on carbon the impact parameter is limited to 
6 < 6.5 fm. We checked that with larger b values the 
evaluated cross sections do not change. In the MC sim- 
ulations one cannot model elastic pion-nucleus reaction. 
The sum over all possible interaction channels gives the 
pion-nucleus reaction cross section. 

Figs HUH] show that the NuWro predictions are in a good 
agreement with the data. 

FT effects can be also used in the secondary interaction 
but their impact on the final results is very small. 



IV. FORMATION TIME/ZONE 

A. Generalities 

The concept of the Formation Time/Formation Zone 
(FT/FZ) was introduced by Landau and Pomeranchuk 
[40| in the context of multiple scattering of electrons pass- 
ing through a layer of material. In the LAB frame the 
FT is given as: 



FIG. 7. Average number of backwards going pions as a func- 
tion of in the NOMAD experiment. 



with a four-momentum pg = [Eg, 0,0, Eq — Mq) and 
the photon by a secondary hadron with a four-momentum 
= {E,pT, ^/E^ -p^~ M2) obtaining: 



t tf 



Eo 



EEo - ^E^-pj,^M^^Ei~M§ 



1 



£2 



1-^ 



(7) 



Pt 



where fj,T is the transverse mass defined = M 

For higher energies E :s> ^j-t, Eq :s> Mq and 

2E 



t 



f 



(8) 



where x = Rantf p^] argued that a further simpli- 
fication a; « is usually well justified and finally in the 
LAB frame: 



2E 



M2 



Pt 



(9) 



and in the hadron's rest frame: 



t = 



E 



■P 



(6) 



where = {E,p) and fc^ = {'-^,k) are four-momenta 
of the electron and the emitted photon respectively, t 
has the interpretation of the minimal time necessary for 
a photon to be created. 

The idea of FT was applied to hadron production by 
Stodolsky [Tsj . who considered the multiproduction of 
mesons by protons passing through a nucleus. In the 
Eq. [5] he replaced the electron by a projectile hadron 



t 



2M 



f,rest 



Pt' 



(10) 



Inspired by this expression Rantf postulated another 
formula for the FT in the hadron rest frame. He kept 
the basic relativistic character of the FT but introduced 
an arbitrary parameter tq to control its size: 



M2 



Trest — Tq 



M2 



Pt' 



(11) 
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The FT defined in Eq. [TT] was implemented in the 
MC event generator DPMJET ^] which later became a 
part of the FLUKA code and was used by the NOMAD 
coUaboration [43|. In the DPMJET cascade model the 
FT is applied to hadrons resulting from all the interac- 
tion modes: QEL, RES and DIS. Following the ideas of 
Bialas [43^ in DPMJET values of FT are sampled from 
the exponential distribution. 

The FT played an important role in the NOMAD anal- 
ysis of the CCQE events ■ They populate mainly one- 
and two- tracks samples. A change of tq modifies the MC 
predictions for the size of both samples: an increase of 
FT makes an impact of the FSI effects on the ejected 
protons smaller and they are more likely to have larger 
momentum with increased probability of being detected 
and populating the two-track sample of events. By ad- 
justing the size of the formation time the values of Ma 
calculated independently from either of the two samples 
of events became almost identical. 

In the above estimations of the FT effect several as- 
sumptions were made which are not necessarily valid at 
lower energies. This is taken into account in the more re- 
cent low energy FLUKA cascade model, called PEANUT 
[44j. In the case of the QEL reaction the concept of co- 
herence length (CL) was proposed to substitute the FT 
effect. 

Derivation of the CL is based on the uncertainty prin- 
ciple arguments: Let be the outgoing nucleon four- 
momentum and q'' — i'-^iQ) the four-momentum trans- 
fer, both in the LAB frame. Because p ■ q is a Lorcntz 
scalar, one can calculate uj ("denotes quantities calculated 
in the nucleon rest frame), the energy transfer in the final 
nucleon rest frame: 



\p-q\ 



\p-q 



— \ujAI\ ^ luj] = 



\p ■ q\ 

M 



From the uncertainty principle uj can be used to esti- 
mate the reaction time in the nucleon's rest frame and 
then in the LAB frame as well. Within that time the 
nucleon is assumed to be unable to re-interact 1451: 



tcL,re 
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(13) 



Surprisingly, the Landau-Pomeranchuk formula Eq. [6] 
is reproduced. 

Among other approaches to give a quantitative eval- 
uation of the FT effects one should mention the SKAT 
parameterization of the LAB frame FZ (46j : 



Is 



KAT 



M 



(14) 



The value of the free parameter was found to be /i^ = 
0.08 ± 0.04 GeV^ based on the experimental data for 
the multiplicity of low momentum (300 MeV/c < p < 



600 MeV/c) protons. This value of fj? agrees also with 
the analysis of the momentum distribution of negatively 
charged mesons in the region p < i GeV/c. 

The SKAT formula can be translated to the following 
value of FT: 
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tsKAT,rest 
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(15) 



Compared to the Rantf formula (Eq. [TT]) the SKAT 
parameterization corresponds to px = but it also in- 
troduces a scale proportional to the hadron mass: tq o 
M//i^. According to the SKAT parameterization FZ is 
identical for pions and nucleons with the same momen- 
tum. At p ^ 1 GeV/c FZ is expected to be ~ 2.5 fm, 
which is of the size of the carbon nucleus. 

In another approach to model FT effects, at a distance 
z from the interaction point one postulates an effective 
(reduced) hadron- nucleon interaction cross section [ivj : 



1 



-zMrn^) /\p\ 



(16) 



with mo « 0.4 GeV. 

Similar description of the FT effects (reduction of the 
cross section) is used in the quantum diffusion model psj : 



'hN ( 



'hN 



< n^kf > 



I'll 



eiin ~~ z) + e{z ~ ih) 



(17) 

where cr^j^ is the free hadron-nucleon cross section, z is 
(12) the distance from the interaction point, k"^ = 0.35 MeV/c 
is the average quark transverse momentum, n = 2, 3 for 
pions and nucleons. The size of FT is determined by It 
which can be evaluated to be 



hi = 2ph 



Af2 



(18) 



where ph and Mh denote hadron momentum and mass, 
and Mn is an intermediate state mass. The precise value 
of AM^ = — is not known and is estimated to be 
between 0.25 and 1.4 GeV^. In ^ the values of 1 GeV^ 
and 0.7 GeV^ are used for protons and pions. 

In the parameterizations given in Eqs ([TBI) and (fT7|) 
smaller effective cross section translates into a larger av- 
erage distance to the first reinteraction point. 

In the case of pions produced via the A excitation and 
decay there is still another natural way to model the FT 
effect. In the INC picture we can treat the A (like in 
the GiBUU approach [s^l) as a real particle propagating 
some distance before it decays. The A lifetime in its rest 
frame is equal jr, with F ss 120 MeV, so in the LAB 
frame one obtains: 
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FIG. 8. Nucleoli Formation Zone in the LAB frame as a 
function of its momentum. 
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FIG. 9. Pion Formation Zone in the LAB frame as a function 
of its momentum 



MT 



(19) 



where i?A is the A energy in the LAB frame. 

We conclude that various approaches lead to similar 
expressions for the FT as far as the dependence on hadron 
momentum is concerned but numerical coefficient and the 
size of the effect can be quite different. 



B. FT models in MC event generators 

Table IIIII summarizes available information about FT 
models in major neutrino MC event generators. 

NEUT uses the SKAT model both for RES and DIS 



MC 


QE 


RES^ 


DIS 


NEUT 




SKAT 


SKAT 


FLUKA 


Coh length 


Rantf 


Rantf 


GENIE 






Rantf-like 


NUANCE 


1 fm 


1 fm 


1 fm 



Note that every MC has its own slightly different definition of 
what does RES and DIS terms mean. 

TABLE III. FT models in MC event generators 



FLUKA uses Eq. [T3]for quasielastic scattering and Eq. 
[Tl]for other processes. 

GENIE uses Eq. [Tl]but in a simplified form, neglecting 



tcc 



To- 



rn 



GENIE assumes the value tq = 0.342/m/c [52]. One 
can check that for the pions the SKAT formula is repro- 
duced. GENIE applies FZ to DIS events and also to the 
non-resonant background events in the RES dynamics 

NUANCE implemented an effective model in which the 
FZ is always equal to 1 fm [sj]. 



C. NuWro FT model 

In NuWro the formation zone effects are implemented: 

• as coherence length (Eq. [T5)) for quasielastic scat- 
terings 

• as A propagation (Eq. [T^ for RES interactions 

• using Ranft model (Eq. Ilip with some parameter 
T for DIS 

There is a smooth transition between the last two mod- 
els at I^ ~ 1.6 GeV. 

With the Ranft model (without the approximation 
Pt — 0) the average value of the FT depends both on 
neutrino energy and on hadron momentum. For a fixed 
value of the neutrino energy, lower hadron momenta typ- 
ically correspond to larger values of the transverse mo- 
mentum and smaller values of the FT. 

In order to fix the value of the parameter t we analyze 
the NOMAD data for the backward moving pions. 



1. Comparison with NOMAD measurement of backward 
moving pions 

To fine tune our model of FZ we make use of the NO- 
MAD experimental data [lH. The average neutrino en- 
ergy in NOMAD is (E^) = 24 GeV and the target com- 
position is dominated by carbon (64.30%) and oxygen 
(22.13%) with small additions of other elements. 
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Scenario 


Without formation zone 


With formation zone 


1 


37.2% 


83.7% 


2 


43.3% 


15.5% 


2a 


22.0% 


8.1% 


2b 


15.6% 


7.4% 


2c 


5.8% 


0% 


3 


2.7% 


0.7% 


3a 


1.9% 


0.6% 


3b 


0.8% 


0.1% 


4 


16.7% 


0.1% 



TABLE IV. Contribution to events with backward n 
different scenarios (description in text) 



from 



We focus on the pion data. Our main observable is 
the average number of backward moving (cos Olab < 0) 
negative pions Btt^ with the momentum p^r between 350 
and 800 MeV/c, as a function of Q^. This observable 
is very sensitive to the FSI effects. Without FSI the 
number of backward moving pions would be very small 
because they appear mainly due to nuclear reinteractions. 
Introduction of the FT makes the FSI effects smaller and 
reduces the number of Btt^. 

Simulations made for various values of t lead us to the 
conclusion that a good agreement with the data is ob- 
tained with r ^ 8 fm/c. Fig. [7] shows average numbers 
of backward moving tt^ reported by NOMAD, and pre- 
dicted by NuWro with and without FZ, as a function of 
Q^. In order to better understand the NuWro perfor- 
mance we analysed various ways in which Btt" appear: 

1. pions are created in the primary vertex and undergo 
quasielastic scatterings during FSI 

2. pions are created during FSI in pion-nucleon inter- 
actions 

(a) single pion production 

(b) double pion production 

(c) triple pion production 

3. pions are created during FSI in nucleon-nucleon in- 
teractions 

4. there are more pion production processes during 
FSI. 

Contributions from the above scenarios to events with 
backward moving tt^ are listed in Table IIVI 

Table rVl shows distributions of the number of Btt^ in a 
single event. The NuWro predictions (with and without 
the FZ) are compared with the Nomad data [ssj . 



2. Comparison with other MC event generators 

Figs [5] and [HI show the values of the FZ in NuWro com- 
pared to other MC neutrino event generators. It is in- 
teresting that in the case of pions various models of the 



< #B7V > 


Data 


NuWro 


Without FT 


With FT 





939617 


921048 


937883 


1 


4238 


22590 


6126 


2 


164 


375 


8 



TABLE V. Contribution to Bvr coming from events with 0, 
1 or 2 pions 
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FIG. 10. K2K: NC Itt" production as a function of tt" mo- 
mentum; the data and NuWro predictions are normalized to 
the same area. 



FZ give very similar results, while in the case of nucleons 
the differences are much larger. 

The NuWro results are given for a specific neutrino 
energy {E^, — IGeV). The FZ grows with E^, due to the 
transverse momentum in the denominator, which goes to 
zero when hadron energy becomes higher. 
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FIG. 11. SciBooNE: NC tt" production as a function 7r° mo- 
mentum; the data and NuWro predictions are normalized to 
the same area. 



10 



2.0 



o 




O 


1 .5 






> 














1.0 


a 








'o 




'dp[ 


0.5 


e 









MB data 








If ^ 




With formation zone 















200 



400 



600 



800 1000 1200 1400 



n momentum [MeV/c] 



FIG. 12. MiniBooNE (neutrino mode): NC ln° production 
as a function tv° momentum. 
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FIG. 13. MiniBooNE (anti- neutrino mode): NC Itt" produc- 
tion as a function tt'' momentum. 



In NuWro the dependence of the FZ on the pion kinetic 
energy is very flat. 



V. APPLICATION: NC Itt" PRODUCTION 
A. Free nucleon NC tt" production 

The data for NC Itt" production cross section on a 
free nucleon target is very scarce. The only such mea- 
surement so far was done in the Gargamelle bubble cham- 
ber. The target was in fact composed of CsHg (90%) and 
CFsBr but the FSI effects were subtracted according to 
the model in (58| . 

In view of large uncertainties in the understanding of 
nuclear effects the results should be treated with some 
caution. We notice also that the data contain a con- 
tribution from the COH reaction. Originally, the results 





Cross section per nucleon (xlO ^^cm^) 


Channel 




NuWro 


free nucleon 


bound nucleon 




0.13 ±0.02 


0.15 


0.12 




0.08 ± 0.02 


0.17 


0.14 


Vfip — >■ Vf^n-K^ 


0.08 ± 0.02 


0.13 


0.11 




0.11 ±0.03 


0.14 


0.13 


u^n fipn'^ 


0.24 ± 0.04 


0.38 


0.36 



TABLE VI. Single pion production cross sections 
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FIG. 14. MiniBooNE (neutrino mode): NC l7r° production 
as a function cos 9. 



were presented as efficiency corrected relative production 
rates in several pion production channels [s^ . The data 
re-analysis was done in [57l |: information about neutrino 
flux was taken into account and the cross sections esti- 
mations were done. 

Table IVII shows the experimental data from [s^l and 
NuWro predictions obtained for neutrinos of energy 
2.2 GeV on free nucleon target and also on nucleons 
bound in ^^C, with Pauli blocking and Fermi motion ef- 
fects taken into account but without FSI effects. We find 
the agreement to be satisfactory. It is also possible to 
compare the NuWro predictions and the data for relative 
contributions from RES and DIS reaction channels. In 
the case of i^p i/pn^ reaction the NuWro predictions 
for the RESiDIS ratio are: 78 : 22 for free and 82 : 18 
for bound nucleons. The experimental data are 80 : 20 
(see Fig. 11 in fs^). 



B. NC Ivr" production on a nucleus 

The recent experimental data for NC 7r° production 
come from three experiments. Basic information about 
them is summarized in Table I VIII 

In the K2K and MiniBooNE (MB) experiments the 
signal was defined as exactly one tt*^ leaving the nucleus 
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FIG. 15. MiniBooNE (anti- neutrino mode): NC ln° produc- 
tion as a function cos 6. 
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FIG. 16. SciBooNE: NC vr" production as a function cos 9. 



target and no other mesons in the final state. In the case 
of SciBooNE (SciB) the signal was defined as at least one 
TT*^ in the final state, with possible other pions as well. 

The experimental signal for Itt*^ production comes 
from: (i) single tt*^ produced at the interaction point in 
the single pion production reaction; (ii) produced in 
double pion production reaction with other pion being 
absorbed; (iii) single tt^ production with charge exchange 
reaction tt^ — >■ tt^ inside nucleus; (iv) primary quasi- 
elastic reaction with being produced due to nucleons 
re-interactions inside nucleus. For the 2Tr'^ production 
the number of possible scenarios is even bigger. 

According to NuWro, most of the Itt" signal events 
(93 — 95%) come from the initial RES single pion pro- 
duction reactions, see Table IVIIII In the case of the MB 
antineutrino flux the contribution is the smallest, because 
the antineutrinos are on average less energetic. Also the 
impact of the FZ is in clear anti-correlation with the av- 
erage fiux energy. 



Experiment 


Beam 


[GeV] 


Target 


Normaliz. 


Measurement 


K2K [59] 




1.30 




relative 


dN/dT„ 


MB [60] 




0.81 


CH2 


absolute 


da / dT^r , 
da /d cos 


MB [60] 




0.66 


CH2 


absolute 


da/dT^ 
da /d cos Ojy 


SciB [61J 




0.81 


CsHs 


relative 


dN/dT^ 
dN 1 d cos Qt, 



TABLE VII. Recent NG tt" production measurements. 



Ghannel 


K2K 


MB V 


MB V 


l7r« ^ Itt" 


93.1% (84.5%) 


93.0% (88.3%) 


94.8% (92.4%) 


no TT — !> Itt'^ 


2.0% (3.2%) 


1.8% (2.4%) 


1.2% (1.6%) 


other TT —5- Itt" 


3.7% (6.8%) 


4.2% (5.8%) 


3.2% (3.9%) 


more tt Itt" 


1.2% (5.5%) 


1.0% (3.5%) 


0.7% (2.1%) 



TABLE VIII. Origin of the events with Itt" in the final state. 
Values in brackets refer to results without FT. 



Table IIXI enumerates what can happen to a tt'' pro- 
duced in the primary vertex due to the FSI effects. One 
can see that pion absorption (the second row) reduces 
the number of NC 7r° events, but the FZ makes the effect 
much smaller. Also the charge exchange reaction (the 
third row) has a significant impact on the final states. It 
is clear that the NC tt*^ production measurement is a very 
good test for the FSI models in MC event generators. 

Table 1x1 shows the composition of the tt" signal in the 
SciB experiment as it is understood by NuWro. The 
second column contains the values reported by the SciB 
collaboration obtained from the MC they used in the data 
analysis (NEUT). 

K2K and SciBooNE did not publish the normalized 
differential cross section. However, flux averaged ratios of 
NCtt" to total CC cross sections were given. In the Table 
IXII we compare both values with the NuWro results. 

Figures [TU] - [121 show the data and NuWro predictions 
for tt'^ momentum distribution in various experiments. 

In the case of the normalized cross section the main 
effect of the introduction of the FZ is the increase of the 
cross section in the pion absorption peak region. The 
effect can be estimated to be 10 — 15%. In the case of 
the K2K measurement the use of the FZ also moves the 



Ghannel 


K2K 


MB V 


MB V 


W - 


-> Itt" 


81.6% (64.0%) 


79.1% (66.9%) 


83.0% (74.5%) 


W - 


no TT 


5.9% (19.3%) 


7.2% (19.2%) 


6.4% (15.9%) 


1^«- 


other TT 


10.1% (11.0%) 


10.2% (10.1%) 


9.6% (7.8%) 


- 


-7> more vr 


2.4% (5.7%) 


2.0% (3.7%) 


1.0% (1.8%) 



TABLE IX. Impact of FSI elects on the events with Itt" in 
the primary interaction. Values in brackets refer to the results 
without FT. 
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Channel SciB MC NuWro (no FT) NuWro (FT) 

Itt" ' 85% 80% 82% 

l7r° + charged tt 11% 16% 14% 

2^ I 4% I 4% I 4% 

TABLE X. Predictions of contribution to tt'' channel from 
SciBooNE MC and NuWro 



NCttVCC 


K2K 


SB 


Data 


0.064±0.008 


0.077 ± 0.010 


NuWro (without FZ) 


0.070 


0.071 


NuWro (with FZ) 


0.079 


0.077 



TABLE XL NCttVCC ratio 



peak of the pion momentum distribution to larger values 
by about 50 MeV/c resulting in much better agreement 
with the data. 

Both MB and SciB experiments provide distributions 
of events versus the cosine of the angle between the neu- 
trino and tt" momenta. Figs [14] - [16] show pions angular 
distributions together with the NuWro predictions. We 
focus on the backward directions because we expect an 
important impact from FZ effects in this kinematical re- 
gion. 

Figs [14] and [15] show that the FZ increases the tt" pro- 
duction in the backward directions but the effect is rather 
small. The reason is that at lower neutrino energies there 
are many backward moving tt'^'s even without FSI effects. 



We checked that only for larger values the FSI become 
the main source of tt^'s and using the FZ reduces their 
number. 

In the case of SciB experiment the NuWro results are 
normalized to the number tt*^ predicted to be in the data. 
In this case using the FZ makes the absolute number of 
backward moving tt^'s little larger, but the effect can 
hardly be seen. 



VI. CONCLUSIONS 

Any comparison to recent NC 7r° production data re- 
quires a computational tool capable of modeling several 
dynamical mechanisms for neutrino-nucleon interaction 
as well as the FSI effects. The NuWro MC event gen- 
erator has all the required physical models implemented 
and we have demonstrated that it reproduces the experi- 
mental results quite well. An important ingredient of the 
NuWro FSI model is the FZ mechanism which even at rel- 
atively small neutrino energies typical for K2K, MB and 
SciB experiments leads to observable effects on the tt'^'s 
in the final state. We hope that our results will be useful 
for better evaluation of the systematic error coming from 
NC 7r° production in neutrino oscillation experiments like 
T2K. 
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